Theory of induced quadrupolar order in tetragonal YbRu2Ge2 
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The tetragonal compound YbRu2Ge2 exhibits a non-magnetic transition at Tb=10.2K and a 
magnetic transition at Ti=6.5K in zero magnetic field. We present a model for this material based 
on a quasi-quartet of Yb 3+ crystalline electric field (CEF) states and discuss its mean field solution. 
Taking into account the broadening of the specific heat jump at To for magnetic field perpendicular 
to [001] and the decrease of To with magnetic field parallel to [001], it is shown that ferro-quadrupole 
order of either Ol or O xy - type are prime candidates for the non-magnetic transition. Considering 
the matrix element of these quadrupole moments, we show that the lower CEF states of the level 
scheme consist of a Ve and a T-? doublet. This leads to induced type of O2 and O xy quadrupolar order 
parameters. The quadrupolar order introduces exchange anisotropy for planar magnetic moments. 
This causes a spin flop transition at low fields perpendicular [001] which explains the observed 
metamagnetism. We also obtain a good explanation for the temperature dependence of magnetic 
susceptibility and specific heat for fields both parallel and perpendicular to the [001] direction. 
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I. INTRODUCTION 

It is well known that some /-electron systems show 
multipole ordering. The phenomenon has attracted much 
attention, since features of multipole order are quite dif- 
ferent from usual magnetic order. As a typical case, 
CeBg shows a kind of antiferro-quadrupole ordering at 
3.4K, whose transition temperature increases with in- 
creasing magnetic field^i 3 -. For transition in NpC>2, an 
octupole ordering is proposed due to experimental results 
of resonant x-ray scattering*^ and NMRZ although a 
cusp at the transition temperature is observed in uniform 
susceptibility^. There are at least two common proper- 
ties between these compounds. At first, these compounds 
have cubic crystal structure. Secondly, the crystalline 
electric field (CEF) ground states of corresponding level 
schemes are quartet states, which are available only in 
cubic systems. It is thought that the quartet state is 
responsible for multipole ordering. 

Recently, some anomalous properties have been ob- 
served in the tetragonal metallic compound YbRu2Ge2 10 . 
In specific heat measurement without magnetic field, 
there are three transition temperatures at To=10.2K, 
Xi=6.5K, and T2=5.7K. It is important that the entropy 
around To obtained by integration of specific heat data is 
very close to Rln4, which means the existence of a quasi- 
quartet state even in the tetragonal system. Applying a 
magnetic field perpendicular to [001] direction, the spe- 
cific heat jump at To broadens, and the peak position 
corresponding to To seems to increase, while T\ and T2 
merge and decrease. Increasing magnetic field further 
above 7T, no anomaly is found. On the other hand, in 
magnetic field parallel to [001], not only T\ and T2 but 
also To decrease with increasing magnetic field. For the 
magnetic susceptibility Xah in magnetic field perpendic- 
ular to [001], no anomaly appears at To, while a cusp is 
observed at T\ for small magnetic fields. Furthermore, 
a metamagnetic transition at higher fields around 2T is 
regarded as spin-flop transition. The magnetic suscepti- 
bility Xc m field parallel to [001] is almost temperature 



independent between To and Ti, and shows flat temper- 
ature dependence below T\ after a slight decrease just 
below T±. Because the value of the paramagnetic effec- 
tive moment 4.5/zb is very close to magnetic moment of 
free Yb 3+ , it is a reasonable assumption that /-hole of 
Yb 3+ is almost localized. Considering the quasi-quartet 
state in a localized picture, some multipole moments will 
be active at each site in the system. From these exper- 
imental data, it has been suggested that To is a kind of 
quadrupole transition, while the phase below T\ is re- 
garded as antiferromagnetic phase with planar staggered 
momentiS. According to Jeevan et al., a change in mag- 
netic structure may happen at T2. We will ignore this 
subtlety in the following and consider only T\ . 

In the theoretical analysis of CeB6, Shiina et al. have 
provided a mean-field approximation for the effective 
Hamiltonian of localized /-electrons belonging to Tg ir- 
reducible representation in Oh point group. In this case 
all multipole moments up to octupole are active*^. The 
relevant multipoles have been classified according to ir- 
reducible representations of the point group in zero mag- 
netic field. Taking into account that some symmetry op- 
erations of Oh point group elements are lost in a magnetic 
field, the multipoles have been reclassified according to 
irreducible representations of the relevant point group in 
the magnetic field. Using these multipoles a mean-field 
approximation has been applied to an effective Hamilto- 
nian to construct the H-T phase diagram. After intro- 
duction of anisotropic interaction of quadrupoles, they 
have obtained a consistent explanation for the anomalous 
ordering in CeBg. It should be noted that this approach 
is promising for multipole ordering not only in CeBg but 
also in TmTe**, where Tm 2+ has the same (4/) 13 elec- 
tronic configuration as Yb 3+ , and NpO?* 3 -. 

In the present work, we apply this approach to inves- 
tigate the phases of YbRu2Ge2. In this system, there 
are some significant differences to CeBg, though a kind 
of quadrupole ordering is expected. At first, the crystal 
structure of YbRu2Ge2 is tetragonal, with point group 
D4h- Second, noting that composition of the quasi- 
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quartet depends on crystalline electric field parameters, 
it is expected that multipolc ordering is also affected by 
the level scheme, while only the size of coupling con- 
stants decides on the favorable multipole ordering in cu- 
bic systems^. Third, considering the present system is 
tetragonal, some multipoles are described only by in- 
duced moments, whose expectation value in the CEF 
ground state vanishe o 14 i 15 . Especially, the second and 
third points bring additional complexity to identify a 
multipole transition. In order to explain the behavior of 
YbRu 2 Ge2, we introduce an effective Hamiltonian in the 
next section. Then, we apply a mean-field approximation 
for the Hamiltonian to identify the non-magnetic order- 
ing state below To- Furthermore, we try to reproduce 
temperature dependences of specific heat and uniform 
susceptibility in magnetic field with anisotropic magnetic 
interaction. Finally, we summarize our results. 



II. EFFECTIVE HAMILTONIAN 

From analysis of uniform susceptibility, effective mo- 
ment is estimated as 4.5/j,b, which is quite close to the 
value 4.54/zb for free Yb 3+ ions. This means that the pic- 
ture of localized hole in the 4/-shell will be reasonable for 
YbR.U2Gc2. First we need to construct CEF level scheme 
of Yb 3+ , to extract relevant multipole moments, and then 
we introduce effective intersite interactions between the 
multipole moments. 



A. CEF term 



addition, the inelastic neutron scattering (INS) gives im- 
portant informations like splitting energy between the 
ground and first excited states in the level scheme. In 
recent INS experiment in YbRu2Ge2, the level scheme is 
reported with the splitting energy of 0.9 meV— . Using 
the splitting energy, we have carried out the fitting of uni- 
form susceptibility. Unfortunately, we could not obtain 
unique CEF level scheme from this procedure. However, 
from reasonable CEF level schemes, we obtain the follow- 
ing common features: (1) The splitting energy between 
the ground and first excited states is about 12K, which 
is estimated by reproducing the entropy obtained from 
specific heat data; (2) the ground and first excited states 
consist of one Tq and one r 7 states, neither two Tq nor 
two r 7 states; and (3) energy splittings of the second 
excited state from the ground state are at least thirty 
times larger than the observed transition temperature 
To=10.2K and the splitting energy between the ground 
and first excited states. Due to the third point, we can 
neglect upper two doublets, if we consider only low tem- 
perature region. The relevant CEF Hamiltonian is then 
given by 

#cef = - -y(fi6nf' l6r i ~ fi7 V fav), (5) 

where Ao is splitting energy from state to state. 
Here, f- is creation operator of /-hole with pseudo-spin 

77 in Kramers-doublet at site i. Although we assume 
that the lower two doublets consist of one r 6 and one r 7 
states, this assumption will be justified during identifica- 
tion of non-magnetic ordered state in YbRu2Ge2. 



The total angular momentum j of Yb 3+ -ion is j = 7/2. 
The multiplet splits into four Kramers-doublets in tetrag- 
onal crystal structure of YbRu2Ge2. In tetragonal point 
group D4h, j — 7/2 multiplet is classified according to 
two-Fg and two-F7 irreducible representations. For two 
doublets belonging to the same r-irreducible represen- 
tation, we call the lower and higher ones rW and r( 2 \ 
respectively, in the following. These states are described 
by linear combination of free ion states \(i)= \ jy) (\n\ < |) 
as follows, 
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where /i is z-componcnt of total angular momentum and 
+ (— ) of left-hand side shows pseudo-spin up (down) in 
Kramers-doublets. 

Usually, CEF parameters are estimated by fitting cal- 
culated uniform susceptibility to the observed one. In 



B. Zeeman term 

The Zeeman term due to the applied field h is given 

by 

Hz = -gjMB X! h ' Ji > ( 6 ) 

i 

where J, gj, and /ib are total angular momentum, Lande 
g-factor of Yb 3+ , and Bohr magneton, respectively. With 
use of second quantization, x- and z-components of the 
angular momentum in r 6 -r 7 subspace are expressed as 

= c 66^66 + c 77^77, CO 
= £66*5(36 + c 77'^77 + c 67 (^67 + Sfe)i (8) 

with a-component of pseudo-spin operator given by 

^tt' — 2 ^ / ^Tri^rm' fr'ri' i (9) 

where o~ a is a-component of Pauli matrix. The coeffi- 
cients c" i are expressed by ct\i and (3\2 as 

cee = 7 - Sa 2 12l c z 77 = -5 + 8^, (10) 
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c& = 4aj 2 , c x 7 = 4 V3/3 12 ^/l-/3? 2 , 

cf 7 = V7y/{1 - a\ 2 ){\ - f3 2 12 ) + %/30a 12 /? 12 . (11) 

Therefore the coefficients ot\2 and /3i2 which determine 
the structure of Tg 1 ' and states are incorporated 
through anisotropic effective ^-factors in Eqs. (I7l8p . We 
comment on the limiting case of ai 2 = 0i2 = 1, which 
are very close to values estimated by fitting of uniform 
susceptibility in III. B.. As we have mentioned above, 
we consider only quasi-quartet consisting of one Tg and 
one IV doublet at each site. As far as the quasi-quartet 
is concerned, multipole moments up to octupole are rele- 
vant as given in Table I. In this sense, the inter-site term 
will be always mapped to j=3/2 quartet system. In par- 
ticular, with ai2 = P12 = 1, the two Kramers doublets 
reduce to | ± |) and | ± |), which are belonging to T 6 
and r 7 irreducible representations in point group. In 
addition the operator J z in j=3/2 quartet system is the 
same as the operator given in Eq. ([7]) with ai 2 — [3\ 2 = 1. 
Therefore, when a magnetic field is applied in [001] direc- 
tion, the present system with a\2 = P12 = 1 is mapped 
to j =3/2 quartet system. However, such mapping is not 
applicable in magnetic field perpendicular along [001], 
due to differences of matrix elements of J x and J y be- 
tween j =3/2 quartet system and the present system with 

"12 = 012 = 1- 

C. Inter-site term 

In the present case, we consider that /-hole localizes at 
each Yb-site. From the simplification mentioned above, 
we have one T 6 doublet and one r 7 doublet. Even in the 
simplification, there are 15 kinds of multipoles at each 
Yb-site. In order to describe the multipoles, we intro- 
duce bases of multipoles <^ belonging to T-irreducible 
representation in point group. In Table. [H we clas- 
sify (fa according to irreducible representations in D4h 
point group, where in addition to S" T ,, we use a /-charge 
operator 

frr-^/Uv (12) 

n 

Since classification of multipoles up to octupole is also 

shown in Table. HI correspondence between multipole 

and </>Jj will be clear. For example, the x-component of 

dipole moment J x is described by a linear combination 

p— 

of 4> n % , which is consistent with Eq. |5]) . 

Now, considering that metallic behavior has been ob- 
served in YbRu2Ge2, effective RKKY interactions be- 
tween the multipoles are present, which are derived 
through the Schrieffer- Wolff transformation from hy- 
bridization term between 4/- and conduction-electrons. 
Noting that the Yb-sitcs in the compound form 
body-centered tetragonal lattice, inter-layer interactions 
should favor ferro-type order, since antiferro-couplings 



would lead to frustration. In the following, we consider 
multipole ordering within a mean-field theory for Yb in 
the body centered tetragonal structure, assuming that 
the ordering takes place either at the zone center q = 
(ferro) or at the zone boundary q=(7r,7r,0) (antiferro). If 
we consider only diagonal term of nearest-neighbour cou- 
plings, the inter-site term of Hamiltonian iJi nt is given by 

= -^EE J n(q)^(-<M(<i)> (1 4 ) 

r,n q 

with (j)n(q) = Yli e_iq Vin> w here j£ is a coupling con- 
stant between and Nq is the number of Yb-sites in 
crystal. In Eq. (14), j£(q) and </>^(q) are Fourier com- 
ponents of coupling constants and 4>\ n , respectively. For 
square lattice, J£(q) is given by 

J£(q) = 2J^(cosq x +cosq y ). (15) 



D. Resulting effective Hamiltonian 

In order to explain low temperature property of 
YbRu2Ge2, effective Hamiltonian used in the following 
is described by 

H c fi = i?CEF + Hz + Hi nt . (16) 

As we have already mentioned, we can determine |A | 
from the entropy. However, in Zeeman term Hz, there 
are two free parameters, which control the weight of 
I ± 1/2) and | ± 3/2) in r 6 1} and states, respectively. 
In addition, we have coupling constants which will 
be estimated in the following sections. Then we apply 
a mean-field approximation for the effective Hamiltonian 
to calculate thermodynamic quantities and the phase di- 
agrams. 



III. ANALYSIS OF TRANSITION AT T 

In this section, we develop a mean-field theory for the 
effective Hamiltonian to analyze non-magnetic transition 
at To. After comparing phase diagrams for all types of 
ferro- and antiferro-quadrupole ordering with experimen- 
tal one, we propose a prefered type of quadrupole order 
in YbRu 2 Ge 2 . 



A. Mean-field approximation for quadrupolar order 

At first, we give mean- field Hamiltonian to determine 
transition line in H-T phase diagram. We assume that 
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TABLE I: Classification of multipoles and relevant bases of 
the multipoles <f>„ according to irreducible representations of 
D4h point group. The first column shows irreducible repre- 
sentation of D4h. The second and third ones describe mul- 
tipoles and cj>n belonging to T irreducible representation, re- 
spectively. Here J, O, and T in the second column are dipole, 
quadrupole, and octupole moments, respectively. The fourth 
column shows expression of corresponding local susceptibility 
of 4>n with splitting energy A from Te to T 7 states. The su- 
perscript ± of irreducible representation expresses the parity 
with respect to time reversal. 
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multipole ordered state is specified by irreducible repre- 
sentation r and ordering wave vector q. From the effec- 
tive model, we obtain easily the mean-field Hamiltonian 

Hmf = Hcef + Hz + H int , (17) 

6 ™ = -j^EE ^n(q)(2<^(q))^(q) 



F,n q 



-(^(q)> ; 



with 



Tr e _ ^ MF / r • • 
Tr e ~H MF /T 



(18) 



(19) 



where T is the temperature. 

We consider three cases, (1) system in zero magnetic 
field, (2) system in magnetic field parallel to [001] direc- 
tion, and (3) system in magnetic field parallel to [100] 
direction, whose point groups are D4h, C^ v , and C2 V , re- 
spectively. In Table (TT] and IIII1 bases of multipoles are 
classified according to irreducible representations of C4 V 
and C2V point group, respectively. In order to develop 



TABLE II: Classification of cf> n according to irreducible rep- 
resentations of C4 V , which is point group of tetragonal system 
in magnetic held parallel to [001]. Here, we note that <jfc are 
bases of multipoles belonging to F irreducible representation 
in D4h point group. The first column shows irreducible repre- 
sentation of C4v The second, third, and fourth ones describe 
quadrupole, c/>JJ , and t/>^ belonging to each irreducible rep- 
resentation in C4 V , respectively. 
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a general formalism, we call bases of multipoles belong- 
ing to T-irreducible representation generically ip^ m an y 
point group G (D4h, C4 V , and Cav)- For D4h point group, 
ipn is equivalent to ^ given in Table. |T] . For non-zero 
field the corresponding lower symmetry point groups C4 V , 
and C2V have basis functions ip!^ that may still be directly 
expressed in terms of the D 4h basis functions as shown 
in Tables II and III. In the following, we discuss only dis- 
ordered, ferro-, and antiferro-ordered states, which is rea- 
sonable since we restrict to nearest-neighbor interaction 
in .Hint. In the disordered state only the fully symmetric 
multipole ^n 1 (®2 m zer0 held) has a non-zero expecta- 
tion value (^^(O)). It leads to the background temper- 
ature dependence of Tq — Tj splitting as shown later. In 
ferro-ordered state belonging to T-irreducible representa- 
tion, (-0^(0)) have finite values in addition to (^^(O)). 
In antiferro-ordered state belonging to T-irreducible rep- 
resentation, allowed expectation values are (^(q)) and 

(Ve(o)>. 

For transition from disordered state to either ferro- 
or antifcrro-multipole ordered state, we usually con- 
sider first- and second-order transitions. In a ferro- (or 
antiferro-) multipole ordered state belonging to T irre- 
ducible representation in point group, expectation values 
of multipole moments (^(0)) (or (-0^(q))) have non-zero 
values. In general, free energy of ordered state is lower 
than that of disordered state below the transition tem- 
perature. The explicit expression of free energy in mean 
field approximation has already been given by Shiina et 
al^. In particular, if the transition is of second-order, 
all order parameters continuously reduce to zero on ap- 
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TABLE III: Classification of (f>„ according to irreducible rep- 
resentations of C 2v , which is point group of tetragonal system 
in magnetic field parallel to [100]. Here, we note that <jfc are 
bases of multipoles belonging to F irreducible representation 
in D4h point group. The first column shows irreducible repre- 
sentation of C 2v . The second, third, and fourth ones describe 



quadrupole, 



and 



belonging to each irreducible rep- 



resentation in C2v, respectively. 
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proaching the transition point. We use linearized mean- 
field equation to determine the second-order transition 
point, by expanding the partition function with respect 
to Hi nt given in Eq. (fTgj) . For transition to ordered state 
specified by irreducible representation T of point group 
and ordering wave vector q it is given by 



(^'(q)) = E^™ 2J ™(q)(^(q)) 



(20) 



with local susceptibility Xnn' defined by 



Ann' 



1/T 



dr (V'in(T)V'iri') -<5r=r 



T 



(21) 



where r is the imaginary time coordinate. The are 
calculated in the limit of (V£(<jO) 0. The linearized 
mean-field equation has non-trivial solution only when 
the maximum eigenvalue of 2;£ Lr J^(q) becomes unity. 
We note that the transition becomes second-order, only if 
the transition point determined by linearized mean-field 
equation (f20)) is the same point where the free energies 
of the two states are equal. 



B. Estimation of coupling constants 

In our mean-field Hamiltonian, there are various pa- 
rameters A , ai2, (3\2 and J^. At first, we consider dis- 
ordered state in zero magnetic field. In this case, only 
allowed multipole is the uniform component of (</> r i ). 



Then, Ao is renormalized as 
1 



A = A 



„ -r -J r ^(0)tanhA 



(22) 



For estimation of Ao and J r i (0)=4J r i , we assume 
A(T = 12K) =12K and A =8K, which gives A(T = 
T = 10. 2K) =12. 8K and J r t =4.3K. By choosing these 
parameter values, we obtain reasonable temperature de- 
pendence of entropy around Tq. In the following calcula- 
tion, we fix these values of A and J r i . 

Now we determine the values of otyi and /3i2- These 
affect the magnetic anisotropy of the uniform susceptibil- 
ity. Experimentally, the uniform susceptibility in mag- 
netic field perpendicular to [001] is quite large compared 
to that in magnetic field parallel to [001]. The magnetic 
anisotropy is successfully reproduced by using the follow- 
ing lower two CEF states: The ground CEF state has al- 
most pure I =F 1/2) character belonging to Te-irreducible 
representation, while the first-excited CEF state is al- 
most pure I ±3/2) belonging to r 7 -irreducible represen- 
tation. Then the lower two CEF states are described by 
a\2 ~ 1 and /3i2 ~ 1, while an ~ 022 ~ OAS. With 
these values of a\2 and /3i2, contribution to J z from the 
first excited T7 state cancels out contribution from the 
ground Tq state, while magnitude of g-factor Cg 6 com- 
ing from the ground Tq state in magnetic field parallel 
to [100] is four times larger than c| 6 . Therefore, this 
CEF level scheme seems to be consistent with the ob- 
served magnetic anisotropy of the uniform susceptibility. 
In the following, we study transitions in the compound 
with these parameter values. 

Furthermore, we have many interaction parameters in 
the inter-site term. Because the highest transition takes 
place at To, each multipole interaction has an upper limit. 
We estimate the upper limit of each coupling constant, by 
assuming that the transition to each ordered state takes 
place at To in zero field. Since we use the linearized 
mean- field equation (|2"0|) , we need to evaluate Xnn'- l n 
point group, the matrix % Lr has diagonal form for 
each r. Each eigenvalue of n-th component in T is given 
by 



A, 1 , = 2j£(q) X 



v.. nn 



(23) 



where expressions for Xnn are summarized in 4th col- 
umn in Table. Q] Specifically, for each of the three 
two-dimensional representations, (J x , J y ), (Tf, T"), and 
(Tf , Tjf) in Table I, we assume that the upper limit of 
coupling constant is independent on n. Together with 
one-dimensional representation, the upper limit of cou- 
pling constant for T irreducible representation J^(q) is 
then given by 



2J c r (q) = 



1 



max x^(T = T ) 



(24) 



In Table [TVl by using the splitting energy A determined 
above, we list the calculated values of l/8max Xnn(T = 



6 



TABLE IV: Coupling constants used for analysis of 
YbRu2Ge2. The first column shows irreducible representation 
of D4h. The second column shows upper limits of magnitudes 
of coupling constants for having rj^ order parameter at To. 
The third column shows coupling constants used in the model 
calculation of Sec. V. Here plus/minus sign for J r denotes 
ferro-/antiferro coupling. The value in parentheses exhibits 

reduced coupling constant between <f> 3 !* . 
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Jo), which is equal to upper limit value of J r with fac- 
tor 1/22=1/8, where z is the square lattice coordination 
number. If we use corresponding upper limit for every 
coupling constant, all transition temperatures to ordered 
states would become degenerate at To. In calculating the 
transition to an ordered state with a primary order pa- 
rameter </>^, we assume that only the coupling constant 
for <^>Jj is at the upper limit. All others are reduced by a 
factor a T , {V ^ T) 

j£(q) = J c r (q), 

jr(q) = ar'J c r '(q) forrvr (25 ) 

where < ar' < 1 is a ratio of actual coupling constant 
J r to the critical value . Present available experimen- 
tal data do not allow a unique determination of coupling 
constants. Only that for the primary ferroquadrupole 
order parameter may be fixed by the transition tempera- 
ture To. In order to obtain a low-field behaviour of Tq(/i) 
which is insensitive to the details of the model we assume 
ar"=0.5 for the remaining coupling strengths. 

C. Order parameter candidates for transition at To 

Before we discuss candidates for transition at To, we 
note a property of lower two CEF states for magnetic 
field parallel to [001]. As we have mentioned in previ- 
ous subsection, the ground CEF state belonging to T 6 is 
lower than the first-excited CEF state belonging to r 7 by 
12. 8K at Tq in zero magnetic field. Since the magnetic 
field parallel to [001] does not break 4- fold symmetry, we 
still can distinguish the and r 7 states. In this case 
their energies are easily obtained as — 4 T gj^B^gh for 
the r 6 Kramers doublet and y =p gj/KBC 77 /i for the r 7 
state. There are two remarkable points in the present 



level scheme; (1) c| 6 has negative sign while c 77 is posi- 
tive, and (2) magnitude of c 77 is three times larger than 
that of Cg 6 . From these facts we can expect a level cross- 
ing in the high field region of the disordered phase, where 
the CEF ground state changes from T 6 in low field to T 7 
in high fields. 

In the following we discuss the transition line of Tq(K) 
in magnetic field parallel to [001]. From experimental 
result, the specific heat data shows that the transition 
temperature of the non-magnetic phase decreases with in- 
creasing magnetic field. As mentioned before, we assume 
that the non-magnetic transition is obtained by ordering 
of quadrupole moment. However, we do not know the 
character of the quadrupole moment among the possible 
0°, Of) O xy , and (O yz , O zx ) cases, and even whether 
the ordering is of ferro- or antiferro-type. In order to 
identify the type of quadrupole ordering, we calculate 
transition lines for all kinds of quadrupole orderings to 
compare with the experimental one. Firstly, since cou- 
pling constant between 0\ is ferro-like as estimated in 
the previous subsection, and ferro-component of 0% is 
allowed even in the disordered state belonging to Y\ in 
C4 V point group (Table ITT|), we exclude any ordering of 
0\ from candidates of the transition at Tq. It does not 
break any symmetry and hence cannot lead to a specific 
heat jump. Secondly, we show transition lines for other 
types of quadrupole orderings in Fig. 1, where for cou- 
pling constants other than the primary order parameter, 
ar"=0.5 is chosen. Considering the behavior of exper- 
imental specific heat, ferro-O^ an( A ferro-O^y orderings 
are consistent with the experiment among all types of 
quadrupole orderings. From calculated results we con- 
clude that transition to ferro-0| (or ferro-O^) ordering 
state is of second-order in low magnetic field region. In 
order to explain a decrease of transition temperature in 
the magnetic field, we show schematic view of magnetic 
field dependence of level scheme with possible transitions 
for 0\ and O xy in Fig. 2. In both 0\ and O xy ordering 
states, the ground CEF state |r 6 +) couples only with 
highest state |r 7 +), therefore, 0\ and O xy order are of 
induced type because the quadrupole expectation value 
in the ground state vanishes. In addition, splitting en- 
ergy between these states increases with increasing mag- 
netic field. Then, since the energy denominator of the 
dominant term of corresponding local susceptibility Eq. 
(|2ip increases, the second-order transition temperature 
decreases with increasing magnetic field. On the other 
hand, for ordering of O yz or O zx , transition temperature 
of these quadrupole orderings increases with increasing 
magnetic field, since splitting energy between the ground 
CEF state and excited state |r 7 — ) decreases. 

We now consider behavior in magnetic field perpendic- 
ular to [001]. From specific heat experiment, it is shown 
that the anomaly of specific heat jump broadens with in- 
creasing magnetic field. This means that the transition 
at T) in zero magnetic field reduces to a crossover in mag- 
netic field perpendicular to [001]. We note that crossover 
does not break any symmetry. Instead of transition tern- 
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perature To, we define a characteristic temperature T* 
as an inflection point of the specific heat divided by tem- 
perature. Let us assume antiferro-quadrupole ordering, 
which breaks at least translational symmetry. Since this 
would cause a specific heat jump, we exclude antiferro- 
quadrupole ordered state from the scenario of crossover. 
On the other hand, ferro-quadrupole ordering in general 
breaks local rotational symmetry, although the transla- 
tional symmetry is preserved. However, ferro-quadrupole 
ordering belonging to I\ -irreducible representation nei- 
ther breaks translational symmetry nor rotational sym- 
metry. For the system in magnetic field parallel to [100], 
0\ and 0\ belong to Ti-irreducible representation, given 
in Table HTT1 We note again that ferro-Ofj ordering does 
not break any symmetry for any field direction and there- 
fore is excluded. Thus, we have one scenario that T* will 
be a crossover temperature from usual high temperature 
region to low temperature region where a considerable 
ferro-0| component appears. 

We now consider magnetic field along [110] direction, 
which is also perpendicular to tetragonal c axis. If we 
repeat similar discussion as for [100] direction, the only 
possible ordering is ferro-O^ belonging to Ti irreducible 
representation in corresponding point group (see Table 
V), instead of ferro-(9|. Unfortunately, we cannot dis- 
tinguish ferro-C>2 horn ferro-O xy ordering, because the 
anisotropy in ab plane is not clarified below To from the 
present experimental data. Therefore, in order to distin- 
guish the type of quadrupole ordering in YbRu2Ge2, de- 
termination of different thermodynamical properties for 
two magnetic field directions, [001] and [110], is required. 
In particlular the possible lattice distortion induced by 
the quadrupole ordering should be investigated. Further- 
more, determination of elastic constant softening would 
be desirable. 

In addition, we consider difference of non-magnetic 
multipole moments between two Tq (two T^) level scheme 
and the Fg-F7 system. In case of Tg-F7 level scheme, non- 
magnetic multipole moments are summarized in Table I, 
II, III, and V, where 0\ and O xy are described by linear 
combination of off-diagonal operators p TT > and S* T , with 
orbital indices r ^ r', respectively. On the other hand, in 
the cases of two and two level schemes, correspond- 
ing orbital off-diagonal operators are not of quadrupole- 
type and belong to Ti and T2 irreducible representations 
in and C^ v point groups, respectively. This is be- 
cause associated Kramers doublets belong to the same 
irreducible representation in such level schemes. Further- 
more, in C2V point group these operators belong to Ti 
and T3 irreducible representations, respectively. Due to 
these differences of irreducible representations from those 
in the Fg-Fy level scheme, no consistent explanation for 
transition at To is obtained for level schemes with two Tq 
or two Ti doublets. Considering these arguments, only 
Tg-r7 level scheme provides a reasonable explanation for 
the transition at To. 

Finally, we summarize this section by proposing candi- 
dates for quadrupolar state below Tq. In this subsection, 




FIG. 1: h-T phase diagrams for quadrupole ordered states. 
The solid line determined by linearized mean-field equation. 
The open circle is given by comparison of free energies be- 
tween disordered state and corresponding quadrupole order- 
ing state. When these instability points are (not) the same, 
the transition is of second-order (first-order). It should be 
noted that h' c (To) = 00 in (a-d), which has been shown in 
Appendix B of Ref. 11. 
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FIG. 2: Schematic view of lower two CEF states in magnetic 
field parallel to [001]. By operators of 0\ and O xy , which 

are equivalent to (f) r 3 ~-^{pm + pro) and 4> r i — ^75 (Sot — 
^75), respectively, the ground Ye state couples only with the 
highest I\ state in low magnetic field region. This leads to 
the induced type quadrupole order below To. The energy 
difference between these two states increases with increasing 
magnetic field. 



we have calculated transition lines of quadrupole ordering 
in magnetic field parallel to [001], and discussed specific 
heat jump in magnetic field perpendicular to [001]. Due 
to these results obtained from different point of views, 
only two appropriate candidates remain. The reasonable 
candidates for non-magnetic transition at To are ferro-O^ 
and ferro-Ozy order parameters. In addition, these order 
are of induced type in the magnetic field parallel to [001]. 
Furthermore, we stress that both ferro-O^ an d ieiro-O xy 
orderings can appear only for the Tg-ry level scheme. 
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TABLE V: Classification of (j)^ according to irreducible rep- 
resentations of C 2v , which is point group of tetragonal system 
in magnetic field parallel to [110]. Here, we note that 0^ are 
bases of multipoles belonging to T irreducible representation 
in D4h point group. The first column shows irreducible repre- 
sentation of C2v The second, third, and fourth ones describe 
quadrupole, , and 0jJ belonging to each irreducible rep- 
resentation in C2v, respectively. 
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IV. DISCUSSION OF THE SECOND 
TRANSITION AT T x 

In the previous section, we have analyzed non- 
magnetic transition at To, and proposed cither fcrro-O^ 
or ferro-Ozj, ordering in system with level scheme. 

In this section, we provide some proposals for the sec- 
ond transition at Ti, which is below To- At first, we 
summarize experimental results below T±. Specific heat 
data show-i£ that a 2 nd order phase transition appears 
at T\. The transition temperature decreases with in- 
creasing magnetic field both parallel and perpendicular 
to [001]. Below Ti, uniform susceptibility data exhibit 
a clear difference between the temperature dependence 
of Xc and Xab, which are the uniform susceptibilities in 
magnetic field parallel and perpendicular to [001], respec- 
tively. Here, % c seems to be independent of temperature 
except for a slight decrease just below Ti, while Xab shows 
a clear cusp at Ti in weak magnetic field. From the dif- 
ference of temperature dependence of %c and Xab it has 
been proposed that below Ti an antiferromagnetic state 
appears with staggered magnetic moment perpendicular 
to [001]. In addition magnetization for field perpendicu- 
lar to [001] has shown a metamagnetic transition. This 
has been regarded as a spin-flop transition in the ab- 
plane. 

From a theoretical point of view, time reversal sym- 
metry must be broken eventually below non-magnetic 
transition temperature To, in order to release remaining 
entropy i?ln2 of the Kramers doublet ground state. In 
this sense, the antiferromagnetic ordering is reasonable. 
However, in addition to dipole (magnetic) moments the 



present model has octupole moments which also break 
time reversal symmetry. Below the quadrupolar transi- 
tion temperature T , the point group reduces from 
to D2h in zero magnetic field, because both ferro-O^ an d 
ferro-O^y ordering break only four-fold rotational sym- 
metry. For dipole and octupole moments, there are four 
one-dimensional irreducible representations in D2h point 
group. Among these four irreducible representations, 
three have respective component of dipole moment J a in 
addition to two components of octupoles. On the other 
hand, the remaining irreducible representation has only 
one octupole component. Let us consider the uniform 
magnetic susceptibility in pure octupole ordered state. In 
very weak magnetic field, it is expected that the uniform 
susceptibility does not considerably decrease below the 
transition temperature of the octupole ordering, because 
there is no dipole order parameter in the state. Since this 
behavior is inconsistent with experimental data of x&b, 
we exclude the pure octupole from candidates for the or- 
dered state below Ti. Therefore, even though we include 
octupole degrees of freedom, the state below Ti seems to 
be inconsistent with pure octupole ordering, but rather 
must be an antiferromagnetic state with a considerable 
magnitude of the dipole moment. 

As candidate below Ti , an antiferromagnetic state will 
be reasonable from the above discussion of the suscepti- 
bility. However, the direction of the staggered moments 
in the ab plane is still unclear. The direction depends 
on whether the ordering quadrupole moment below To is 
0\ or O xy because their remaining symmetry axis is dif- 
ferent. Let us consider ferro-0| order in magnetic field 
parallel to [100]. According to Tablc lLTTl thc allowed direc- 
tions of magnetic moments in the system are J x e x , J v e y , 
and J z e z , where e x , e y , and e z are unit vectors parallel 
to [100], [010], and [001] directions, respectively. Further- 
more the direction of staggered moments in the antifer- 
romagnetic state is perpendicular to [001] axis according 
to the experimental uniform susceptibility. Therefore, for 
magnetic transition within the 0\ phase, the dipole order 
parameter should be described by J x e x or J v e y . Like- 
wise within the O xy phase in magnetic field parallel to 
[110], it should be J x e x ± J v e y . This means that deter- 
mination of staggered moment direction can distinguish 
between quadrupole order of ferro-0| or ferro-O^ type. 

Now we consider the metamagnetic transition, which 
has been observed in magnetic field perpendicular to 
[001], starting from an antiferromagnetic state with pla- 
nar staggered moment. In the present scenario, the an- 
tiferromagnetism is regarded to appear below ferro-Of 
or ferro-O^y ordering temperature. In these cases the 
point group reduces from D4h to D2h in zero magnetic 
field. Consequently two-dimensional irreducible repre- 
sentation in D4h, which involves two planar components 
of magnetic moment, reduce to the direct sum of two 
one-dimensional irreducible representations, where each 
has one planar component of magnetic moment. 

Then it is expected that exchange coupling constants 
between planar components of magnetic moments de- 
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pend on the in-plane direction. These effective cou- 
pling constants are due to RKKY mechanism, therefore 
their anisotropy will be induced by reconstruction of con- 
duction electron states in the ferro-quadrupolar ordered 
phase. This exchange anisotropy induced by quadrupole 
order is the origin of the spin flop transition for H _L[001]. 



V. MEAN-FIELD ANALYSIS OF BOTH 
ORDERED PHASES 

From previous discussions, we have two candidates for 
successive transitions; one scenario is given by ferro-O^ 
ordering for the first transition at To before the second 
transition at T± to antiferromagnetism with staggered 
moment J x e x or J v e yi and another is ferro-O^ ordering 
for transition at To before transition at T\ to antiferro- 
magnetism with staggered moment J x e x ± J v e y . In the 
following, we assume ferro-O^ ordering for the transition 
at To, because we cannot distinguish these two possibili- 
ties from available experimental data. 



A. Re-estimation of coupling constants 

In the previous section, we have estimated parame- 
ter values of model Hamiltonian, such that magnetic 
anisotropy of uniform susceptibility and non-magnetic 
transition temperature are reasonably reproduced. How- 
ever, in that stage, we have assumed that all signs of 
coupling constants are the same. In the present case, we 
are considering that the system shows ferro-quadrupole 
ordering before antiferromagnetic transition. Therefore, 
we should estimate coupling constants with assumption 
of ferro-quadrupole transition at To and antiferromag- 
netic transition at T\. Among these coupling constants, 
J r i and J r 3 are not changed from values used for Fig. 
1(a). On the other hand, we choose value of J r s as 
antiferromagnetic transition appears at Ti=6.5K in zero 
magnetic field. In addition, magnitudes of other coupling 
constants are chosen to be small with a r =0.25 (r'=r^~, 
Ig , rr, and T7), so that unobserved phases are 
suppressed which would appear if a r were larger. Fur- 
thermore, considering that coupling constant between x- 
components of magnetic moment is different from corre- 
spondence between y-components of magnetic moment in 
ferro-O^ ordering state, we will reduce magnitude of cou- 

pling constant between <p 3 y . In Table HVl we summarize 
the revised values of coupling constants. 



B. Phase diagram 

At first, we construct H-T phase diagram in magnetic 
field parallel to [001]. As we have mentioned in §111. A, 
in order to find the transition line, we use two kinds of 
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FIG. 3: h-T phase diagram in magnetic field parallel to [001]. 
The solid line shows level crossing temperature. The dashed 
and dotted lines are instability lines of disordered state to 
ferro-O^ state and antiferromagnetic state with staggered mo- 
ment parallel to [100], respectively. The open circle corre- 
sponds to transition point from disordered state to ferro-O^ 
ordered state. The solid circle corresponds to transition point 
from paramagnetic state (including ferro-O^ ordered state) to 
coexistent state of both ferro-O^ moment and staggered mag- 
netic moment parallel to [100]. These circles are determined 
by comparison of free energies. 



procedures; one is given by linearized mean-field equa- 
tion, while the other is determined by comparison of free 
energies of different states. In Fig. 3, we show calcu- 
lated H-T phase diagram in magnetic field parallel to 
[001]. In the high field region of the phase diagram, level 
crossing between and states is obtained in terms 
of larger magnitude and negative sign of (/-factor of the 
T7 state. The transition at To is of second-order, and 
To decreases with increasing magnetic field. The transi- 
tion at Ti is second-order transition from ferro-O^ phase 
to coexistent phase of ferro-0| and antiferromagnetism. 
Here, dashed line shows instability of paramagnetic state 
to ferro-0| state, while dotted line is instability line of 
paramagnetic state to antiferromagnetic state with stag- 
gered moment parallel to [100]. The dashed line is very 
different from instability line shown in Fig. 1, since cou- 
pling constant between field induced is changed from 
ferro-coupling to antiferro-coupling. Comparing these in- 
stability lines with the transition line to the coexistent 
phase, ferro-O^ and antiferromagnetism are cooperative 
to each other and stabilize the coexistent phase. 

In Fig. 4, we show calculated H-T phase diagram in 
magnetic field parallel to [100]. In this figure, To is not 
a transition temperature but the crossover temperature 
from usual paramagnetic phase in high temperature re- 
gion to disordered phase with considerable ferro-O^ mo ~ 
ment in low temperature region, where the crossover tem- 
perature is determined by inflection point of temperature 
dependence of C'/T. The crossover temperature increases 
with increasing magnetic field. In order to consider low 
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FIG. 4: h-T phase diagram in magnetic field parallel to [100]. 
The open circle shows crossover temperature determined by 
inflection point of temperature dependence of C/T. The dot- 
ted and dashed lines are instability lines of disordered state 
to antiferromagnetic state with staggered moment parallel to 
[100] and parallel to [010], respectively. The up-triangle and 
down-triangle correspond to transition point from disordered 
state to antiferromagnetic state with staggered moment par- 
allel to [100] and parallel to [010], respectively, which are de- 
termined by comparison of free energies. Inset shows magne- 
tization in field parallel to [100]. 



temperature state, we take into account reduced coupling 
constant between J y , which is mentioned in §.IV. By 
the anisotropic magnetic interaction, we have two an- 
tiferromagnetic phases. Here, low field phase has stag- 
gered moment parallel to the magnetic field, while high 
field phase has perpendicular component to the magnetic 
field. Therefore, spin-flop transition is obtained, where 
the transition between these phases is of first-order. If 
we use the same coupling constant as the one between 
J x , low field antiferromagnetic phase disappears. In the 
inset of Fig. 4, metamagnetic transition corresponding 
to the spin-flop is obtained. 



pression of internal energy E\, as follows: 

£l = EE £ »<C> "EE ^-jn<^n><4«>> ( 26 ) 
j n i#j r,n 

where the first term of right hand side exhibits contribu- 
tions from CEF term and Zeeman term. Then, specific 
heat is given by temperature derivative of the internal 
energy as 



^ = EE* 



1 dT 



-EE^-in(C)^, (27) 

m r.n 

with equation of temperature derivative of (V>jn) 



dT 



= V — Y ri ,' 



i(?y)r',n' 



T dT 'X-i n ' n '\ ' 



where xlrJn * s g iven by the sublattice dependent expres- 



sion 



Am'n 



1/T 



dr (<,(r)0 



T 



(29) 



C. Specific heat 



Using these expressions, we calculate the specific heat 
of the system in a magnetic field. In Fig. 5(a) and 5(b), 
we show temperature dependences of specific heat for 
field parallel to [001] and [100] directions, respectively. 
In Fig. 5(a), there are two specific heat jumps corre- 
sponding to transitions to ferro-0| state at To and to 
coexistent state at T\, In addition, field dependence of 
the specific heat is very weak in low field region. In Fig. 
5(b), only specific heat jump is due to transition to anti- 
ferromagnetic state at T%, while hump structure related 
to the crossover temperature at To broadens with increas- 
ing magnetic field. This means that the ferro-0| moment 
is induced by magnetic field parallel to [100] as shown in 
Table Hill and no symmetry breaking at To takes place for 
H > 0. 



In specific heat data, there are some anomalous fea- 
tures in the temperature and field dependences. In mag- 
netic field parallel to [001], specific heat jumps are ob- 
served at To and Ti, and these transition temperatures 
decrease with increasing magnetic field. In magnetic field 
perpendicular to [001], the anomaly at the non-magnetic 
transition reduces to a hump and broadens with increas- 
ing magnetic field. In order to compare with the experi- 
mental data, we calculate specific heat in magnetic field, 
based on the mean-field solution. 

In order to calculate specific heat, we first provide ex- 



D. Uniform susceptibility 

The experimental data of uniform susceptibility ex- 
hibits characteristic properties. In the magnetic field per- 
pendicular to [001], uniform susceptibility does not show 
anomaly at To, except for the slight increase of mag- 
nitude of the temperature derivative at To in low field 
region. On the other hand, temperature dependence of 
uniform susceptibility in magnetic field parallel to [001] 
has a plateau- like behavior between To and T\. In ad- 
dition, the temperature dependence is almost insensitive 
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FIG. 5: Temperature dependences of specific heat in magnetic 
field parallel to [001] (a) and parallel to [100] (b). The solid 
and dotted lines correspond to IT and 4T, respectively. 
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FIG. 6: Temperature dependences of uniform susceptibilities 
in magnetic field parallel to [001] (a) and parallel to [100]. The 
solid and dotted lines correspond to IT and 4T, respectively. 



to the magnetic field up to 3T. In order to analyze uni- 
form susceptibility, we give expressions of the quantity in 
magnetic field parallel to [001] and [100]. For direct com- 
parison with experimental data in finite magnetic field, 
magnetization divided by magnitude of the field is used 
as uniform susceptibility, instead of the Kubo formula 
of susceptibility from linear response theory. Then, uni- 
form susceptibility Xc in magnetic field parallel to [001] 
is given by 



Xc 



5.7 MB 



(r 



(30) 



while uniform susceptibility in magnetic field parallel to 
[100] is similarly obtained from 



Xa = — — (J >, 



(31) 



where J z and J x are given by eqs. (6) and (7), respec- 
tively. 

In Fig. 6(a), we show calculated temperature depen- 
dence of uniform susceptibility in magnetic field parallel 
to [001]. Anomalies due to two transitions at T and T\ 



are obtained. Here, we note that plateau in tempera- 
ture dependence between To and T\ is due to moderate 
coupling constant between octupoles cf) T i , which are in- 
duced by the magnetic field, while temperature indepen- 
dent behavior well below T\ reflects existence of stag- 
gered moment perpendicular to the magnetic field. In 
addition, the uniform susceptibility hardly has field de- 
pendence at least in the low field region. In Fig. 6(b), 
we show temperature dependence of uniform susceptibil- 
ity in magnetic field parallel to [100]. From the figure it 
shows slight upturn around crossover temperature T . In 
magnetic field of IT, it shows a cusp at transition tem- 
perature to antiferromagnetism with staggered moment 
parallel to the magnetic field. With increase of magnetic 
field to 4T, the uniform susceptibility shows tempera- 
ture independent behavior below transition temperature 
to antiferromagnetism with staggered moment perpen- 
dicular to the magnetic field. The different of behavior 
below antifcrromagnetic transition temperature is due to 
the spin-flop transition, as shown in Fig. 4. 
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VI. DISCUSSION AND SUMMARY 

Before we summarize, we would like to comment on a 
few points. In the previous section, we have calculated 
specific heat, uniform susceptibility, and phase diagram 
with assumption of ferro-0| ordering at To. Compar- 
ing our results with experimental data of YbRu 2 Ge2, the 
present model not only explains the experimental data 
qualitatively, but also gives quantitative agreement in 
specific heat jumps. Therefore, ferro-O^ ordering at To 
above antiferromagnetic transition temperature T\ will 
be one of promising candidate for non-magnetic transi- 
tion of YbRu2Ge2- However, we have proposed either 
ferro-O^ or iervo-O xy ordering for the transition at To 
in §. III. Considering that both candidates are ferro- 
ordering states of quadrupoles, in order to identify the 
non-magnetic state among the two candidates, it will be 
useful to carry out ultrasonic and x-ray scattering exper- 
iments , because ferro-quadrupole couples with lattice 
distortion. From this point of view, it is desirable to 
confirm the crystal structure below To. 

Related to property of non-magnetic ordering state, it 
has been recently reported that Ru-NQR spectrum may 
not be affected by the non-magnetic transition^. With 
respect to this result, taking into account positions of 
Yb and Ru ions, if dominant quadrupole of Ru-nucleus 
is either OS? or 0\ > ferro-O^j, ordering of Yb 3+ does not 
change Ru-NQR spectrum by the transition. On the 
other hand, if quadrupole of Ru-nucleus is O xy , ferro-Of 
ordering does not change the spectrum by the transition. 
Furthermore, if quadrupole of Ru-nucleus is either O yz 
or O zx , both ferro-0| ordering and ferro-O^j, ordering 
do not change the spectrum. In order to identify the 
type of quadrupole ordering, it is required to clarify the 
quadrupole of Ru-nucleus. 

Furthermore, in recent uniform susceptibility data, it 
is shown that behavior of uniform susceptibility in mag- 
netic field parallel to [100] is similar as that in magnetic 
field parallel to [110], and these uniform susceptibilities 
have finite value at very low temperature. Based on 
the present model with assumption of either ferro-O^ or 
ferro-Ozy ordering at T , one of these uniform suscepti- 
bility will vanish at OK in small magnetic field region by 
development of staggered magnetic moment parallel to 
the magnetic field. Considering that finite values are ob- 
served for both uniform susceptibilities in very low tem- 



perature region, the ordering wave vector of magnetic 
moments will be an incommensurate one. Therefore, it 
is desirable to carry out neutron scattering experiment 
to clarify the magnetic state of the compound. 

Finally, we comment on effect of multipolar fluctu- 
ations. The effect has been investigated in CeBg by 
Shiina— . In the paper, it has been shown that the mul- 
tipolar fluctuation are enhanced by approaching the sys- 
tem to the SU(4) symmetric limit. In YbRu2Ge2, it is 
considered that the tetragonal anisotropy like splitting 
energy A between T§ and T7 states breaks the SU(4) 
symmetry inherently. Therefore we do not expect that 
the multipolar fluctuations change qualitatively behav- 
iors suggested by the mean-field theory. 

In summary, in order to explain properties of 
YbRu2Ge2, we have introduced a quasi-degenerate lo- 
calized model consisting of CEF term, Zeeman term, 
and exchange term of multipoles. Classifying multipolcs 
according to irreducible representations of correspond- 
ing point group, we have developed a mean-field theory 
for the model. Considering that the specific heat jump 
broadens with increasing magnetic field perpendicular to 
[001], we have proposed that for the non-magnetic tran- 
sition, only ferro-0| and ferro-O^y orderings are possible 
candidates. Furthermore, it has been shown that these 
ferro-quadrupole orderings are only available and essen- 
tially of the induced type, when the lower two CEF states 
consist of one Fq and one doublets in zero magnetic 
field. With assumption of ferro-0| ordering at To, we 
have calculated specific heat, uniform susceptibility, and 
phase diagram, where anisotropic exchange interaction 
between planar components of magnetic moments is in- 
troduced as an effect of the ferro-quadrupole ordering. 
The calculated results have been shown to explain exper- 
imental data consistently. In order to clarify the property 
of YbRu2Ge2 completely and refine the set of coupling 
constants, it is desirable to carry out more detailed exper- 
iments like elastic constant measurements and neutron 
diffraction in applied magnetic field. 
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